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We characterize the reflectance peak near the Brewster angle for both an interface between two 
dielectric media and a single slab. To approach this problem analytically, we approximate the 
reflectance by a first-order diagonal Pade. In this way, we calculate the width and the skewness of 

> . 

the peak and we show that, although they present a well-resolved maximum, they are otherwise 

o 

not so markedly dependent on the refractive index. 
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I. INTRODUCTION 



As it is well known, the behavior of the light reflected at the boundary between two 
dielectrics drastically depends on whether the associated electric field lies parallel (p polar- 
ization) or perpendicular (s polarization) to the plane of incidence. This is quantified by 

fl 

the Fresnel formulas which also confirm that at a certain angle 9 b the p-polarized field 
drops to zero. This angle is commonly referred to as the Brewster angle and, in spite of its 



simplicity, it is a crucial concept in the optics of reflection 2|. 



In fact, this notion is at the heart of a number of methods for optical thin-film charac 



terization. These include Brewster-angle microscopy 3, y, |5j , which furnishes a direct visu- 
alization of the morphology of the film, or the Abeles Brewster-angle method 



a direct visu- 



for determining the refractive index. They are quite popular because are noncontact, easy 
to set up and simple to use. The basic idea can be concisely stated as follows: coated and 
bare regions of a substrate are illuminated with p-polarized light and the angle of incidence 
is scanned for a reflectance-match angle. Visual observation is often used because it reveals 
local variations in film thickness and in homogeneity of refractive index. 

These techniques rely on the dependence of reflectance near the Brewster angle on the re- 
fractive indices involved. Of course, one could rightly argue that all the relevant information 
is contained in the Fresnel formulas. However, given their nonlinear features, no physical 
insights on the behavior around 9 b can be easily inferred, nor an analysis of the error sources 



can be easily carried out, except by numerical methods 



111 ] . Some qualitative comments can 



be found scattered in the literature: for example, it is sometimes argued that the amplitude 
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reflection coefficient for p polarization is approximately linear around 9b |12i |. Nevertheless, 
we think that a comprehensive study of the reflectance near the Brewster angle behavior is 
missing and it is precisely the main goal of this paper. 

To this end, we propose to look formally at the reflectance as a probability distribution 
and focus on its central moments, which leads us to introduce very natural measures of the 

n 

width and the skewness of that peak through the second and third moments [l3j. Since 
all these calculations must be performed only in a neighborhood of 9 b, we replace the 



exact reflectance by a Pade approximant 14 1: apart from the elegance of this approach, 
it is sometimes mysterious how well this can work. In addition, we can then compute the 
relevant parameters in a closed way and deduce their variation with the refractive indices. 
This program is carried out in Section |H] for a single interface and in Section 11111 for a 
homogeneous slab. Finally, our conclusions are summarized in Section HVI 

II. BEHAVIOR OF THE BREWSTER ANGLE AT AN INTERFACE 

Let two homogeneous isotropic semi-infinite dielectric media, described by real refractive 
indices uq and n%, be separated by a plane boundary. We assume an incident monochromatic, 
linearly polarized plane wave from medium 0, which makes an angle 9$ with the normal to 
the interface and has amplitude This wave splits into a reflected wave in medium 
0, and a transmitted wave in medium 1 that makes an angle 9\ with the normal. The 
angles of incidence 9q and refraction 9\ are related by Snell's law. 

The wave vectors of all waves lie in the plane of incidence. When the incident fields are p or 
s polarized, all plane waves excited by the incident ones have the same polarization, so both 



basic polarizations can be treated separately. By demanding that the tangential components 
of E and H should be continuous across the boundary, and assuming nonmagnetic media, 
the reflection and transmission amplitudes are given by 1| 



(r) 

Es n cos #0 — ni cos 6\ 
n cos #0 + ni cos 6\ 



(2.1a) 
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(2.1b) 



Ep rii cos 6 — n o cos 6\ 



£jW m cos 8 + n cos 9\ 
These are the famous Fresnel formulas, represented in Fig. 1 and whose physical content is 
discussed in any optics textbook. The Brewster angle occurs when r p = 0, which immediately 
gives the condition 

6 B = tan" 1 n, (2.2) 

where n = ni/n is the relative refractive index. Without loss of generality, in the rest of 
this section we assume Uq = 1 and then n — n\. Also, we deal exclusively with p polarization 
and drop the corresponding subscript everywhere. 
To treat the reflectance 

m=\r\ 2 (2.3) 
near the Brewster angle it proves convenient to use 

$ = 6 -6 B , (2.4) 

which is the angle of incidence centered at 6b- Since we are interested only in a local study, 
we take into account only a small interval around 6b', that is, angles from — Ai9 to +Ai9. 
The length of this interval is largely an arbitrary matter: we shall take henceforth A$ =15°, 
although the analysis is largely independent of this choice. 
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FIG. 1: Plot of the amplitude reflection coefficients r p and r s as functions of the angle of incidence 
#o (in degrees) for an interface between air (no = 1) and a homogeneous, isotropic dielectric medium 
(m = 1.75). 

To quantify the peak near the Brewster angle, we treat the reflectance as a probability 



distribution in the interval [— A$, +A^] and borrow some well-established concepts from 



statistics 



13| . In this manner we define 



(2.5) 



For a symmetric peak i? = 0, while the fact that $ ^ reveals an intrinsic asymmetry. 



The central moments are 

^ = r+A^ > ( 2 - 6 ) 

and, as it happens for they are functions of the refractive index n. The second moment 
is a measure of the width of the distribution, while the third moment can be immediately 
related with a lack of symmetry. More concretely, we take the width and the skewness of 
the Brewster peak as 

w = 2^, (2.7a) 



7 = (2.7b) 



Of course, these parameters can be computed numerically by using the Fresnel formulas. 
However, since we are working in a small neighborhood of 9b, it seems more convenient 
to work with some local approximation to r (and therefore to 9t). Here, we resort to Pade 
approximants (to maintain the paper as self-contained as possible, in the Appendix we briefly 
recall the essentials of such an approach). After some calculations, the approximations to 
the reflection coefficient r turn out to be 

rim = (igv (28a) 



[2/0] = (! - n ^ 
2n 3 

' " "" "' --' y2 (2 8b) 



4n 6 

2n 3 - (n 4 +^ + 2)i?' 

whence one immediately gets the corresponding reflectances as 9^ L//M 1 = |r[ L//M l| 2 , which 
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FIG. 2: Reflectance (for p polarization) of the same interface as in Fig. 1, but calculated in terms 
of the local angle of incidence $ (in degrees) for an interval ±A$ = 15°. The exact reflectance 
is represented by the thick continuous line, while the broken lines are the Pade approximants in 
Eqs. 1)2.8(1 . The corresponding orders [L/M] are labeled in the right axis. 

have been plotted in Fig. 2. As we can see, the term 9^ 1//0 ' reproduces the exact reflectance 
only in a very small interval. The next polynomial approximation 9^ 2 /^ improves a little bit 
the situation, but fails again as soon we are, say ±10° away from 9b- On the contrary, the 
diagonal Pade IH' 1 '' 1 ' fits remarkably well with the exact behavior. We thus conclude that 
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FIG. 3: Width (in degrees) of the reflectance peak near the Brewster angle as a function of the 
refractive index n for the same interface as in Fig. 1 (lower curve) and for a slab of index n and 
thickness d = A/4 (upper curve) . 

fftl 1 / 1 ] provides an excellent approximation and we use it in subsequent calculations. 

An additional and remarkable advantage of D^ 1 / 1 ' is that the central moments can be 
expressed in a closed analytical form. In Fig. 3 we have plotted the width w as a function of 
the refractive index n (lower curve). We have also computed w by numerically integrating 
the Fresnel formulas: no appreciable differences can be noticed. The width has a maximum 
that can be calculated by imposing dw/dn = 0, which immediately gives n = Beyond 
this value, w decreases almost linearly with n. However, for the range of indices plotted in 
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FIG. 4: Skewness 7 of the reflectance peak near the Brewster angle as a function of the refractive 
index, for the same cases as in Fig. 3. The lower curve also corresponds to the interface, while the 
upper one if for the slab. 

the figure, this variation is smooth, with a total change in w of around 3°. 

In Fig. 4 we have plotted the skewness 7 in terms of n. This parameter is always negative 
what, for a peak, means that the left tail is less pronounced than the right tail. There is 
again a maximum in the skewness that can be calculated by imposing d'y/dn = 0. However, 
one can check that ~ 1 for all the values of n, so that 7 ~ —8/w 3 : this maximum 
coincides then with that of w. Apart from a scale factor, 7 shows an aspect quite similar to 
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that of w. 



III. BEHAVIOR OF THE BREWSTER ANGLE AT A SINGLE SLAB 

We focus our attention now on a homogeneous, isotropic dielectric slab of refractive index 
n and thickness d imbedded in air and illuminated with monochromatic light (the case when 
the slab lies on a substrate can be dealt with much in the same way). A standard calculation 



gives the following amplitude reflection coefficient for the slab 1£| 

r[l-exp(-»2/?)] 
l_ r 2 exp (_, 2 /3)' 

where r is the Fresnel reflection coefficient at the interface air-medium (which again will be 
considered for p polarization only) and (3 is the slab phase thickness 

(3 = — d^n 2 - sin 2 O . (3.2) 
A 

Here d is the thickness, A the wavelength in vacuo of the incident radiation and 9q the angle 
of incidence. This coefficient R presents a typical periodic variation with the slab phase 
thickness (3. Apart from that dependence, it is obvious that R — when r = 0; i. e., 
precisely at the Brewster angle for the interface air-medium. However, the form of this 
reflectance peak, calculated now as 

D\=\R\ 2 , (3.3) 

differs from the case of a single interface, since ()3.1|) is more involved than (|2.1|) . Although 
the analysis can be carried out for any value of d such that R ^ 0, for definiteness and to 
avoid as much as possible any spurious detail, we take d — X/A. In Fig. 5 we have plotted the 
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FIG. 5: Reflectance (for p polarization) of a slab of thickness d = A/4 in terms of the local angle 
of incidence for various values of the refractive index (labeled in the right axis). 

exact form of the reflectance for this thickness, as a function of the local angle of incidence, 
defined as in Eq. (|2.4jl . for several values of n 

Instead of Eq. 1)3.1)1 . we use again a Pade approximant of order [1/1], which can be written 



as 



R [1/1] = (3-4) 
1 + qtf y 1 
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where 



Pi = i 




sin u, 



(3.5a) 



v 



n 4 + n 2 + 2 



(3.5b) 



9i = 



n 3 e w sin z/ 



2n 3 



and 



27m 2 d 



(3.6) 



With this approximation, that still reproduces pretty well the exact reflectance, we have 
analytically obtained w and 7 for the slab. The results are plotted in the upper curves of 
Figs. 3 and 4, respectively. For high values of n, the reflectance peak for the slab essentially 
coincides with that for the interface. Differences are only noticeable for small values of n. In 
fact, the maximum of w can be evaluated by imposing dw/dn = 0, whose numerical solution 
gives n ~ 1.45. The behavior of 7 is quite similar, because here we have also that /x 3 ~ — 1 
for all the values of n. 

IV. CONCLUSIONS 

In summary, we have presented a simple and comprehensive treatment of the reflectance 
near the Brewster angle. By combining the notions of width and skewness with Pade ap- 
proximants, we have fully characterized this reflectance peak. We hope that these results 
fill a gap for a full understanding of the Brewster angle. 
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APPENDIX: PADE APPROXIMANTS 

A Pade approximant is a rational function (of a specified order) such that its power series 
expansion agrees with a given power series to the highest possible order. In other words, let 
us define 



The rational function B\ L ^ M \x) is said to be a Pade approximant to the function f(x), which 
has a Taylor expansion at x = 




(Al) 



where Pl(x) and Qm{%) are polynomials of degrees L and M, respectively: 



P L (x) = p + pix + . . . + Plx l , 



(A2a) 




(A2b) 




(A3) 



fc=0 



if 



(A4a) 




(A4b) 
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for k = 1, 2, . . . , L + M. These conditions furnish L + M + 1 equations for p > • • • ,Pl and 
gi, . . . , 5m- The coefficients can be found by noticing that Eqs. (jA4|) are equivalent to 

/(*) - = °> (A5) 

up to terms of order L + M + 1. This gives directly the following set of equations: 

/o = Pa, 

fi + Mi =Pi, 

/a + + /o<?2 = P2, 



/l + fi-iQi + • • • + /l-m+iQm — 0, 



/i + + • • • + h-M+lQM = 0, 

gi+Af + fL+M-XQl + ■ ■ ■ + /l<?m = 0. (A6) 

Solving this directly gives 
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R [L/M] { 



X 



J'l-m+i 


f L-M+2 




Il+i 


h 
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(A7) 



Nevertheless, experience shows that Eqs. f|A6|) are frequently close to singular, so that one 







should solve them by e. g. a full pivotal lower-upper triangular (LU) decomposition 

By contrast with techniques like Chebyshev approximation or economization of power 
series, that only condense the information that you already know about a function, Pade 



lfj. We 



approximants can give you genuinely new information about your function values 
conclude by noting that, for a fixed value of L + M, the error is usually smallest when L = M 
or when L — M + 1. 
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